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Abstract 

We consider the double-radiative decays of heavy-light QED and 



' QCD atoms, fi + e~ — > 77 and B® — > 77. Especially, we take under 

■ scrutiny contributions coming from operators that vanish on the free- 

quark mass shell. We show that by field redefinitions these operators 
are converted into contact terms attached to the bound state dynam- 
5_i ■ ics. A net off-shell contribution is suppressed with respect to the effect 



of the well known flavour-changing magnetic-moment operator by the 
bound-state binding factor. The negligible off-shellness of the weakly 
bound QED atoms becomes more relevant for strongly bound QCD 
atoms. We analyze this off-shellness in model-approaches to QCD, 
one of them enabling us to keep close contact to the related effect 
in QED. We also comment on the off-shell effect in the correspond- 
ing process Bd - > an d discuss possible hindering of the claimed 
beyond-standard-model discovery in this decay mode. 

PACS: 14.60.Ef, 14.40.Nd, 12.15.Mn, 12.15.-y, 13.10. +q 
Keywords: bound state, off-shellness, flavour- changing transition, rare 
decays 



1 Dedicated to Prof. Ksenofont Ilakovac on the occasion of his seventieth birthday 

2 j. o.eeg@fys.uio.no 

3 kkumer@phy.hr 

4 picek@phy.hr 



1 Uses of a comparative study of the off-shell 
effects in QED and QCD atoms 

Off-shell effects are known to be quite elusive. The most famous measured 
effect, dubbed Lamb shift [||,|2[], appears in atomic physics. It is represented 
by the atomic level shift on account of a tiny difference in the self-energies 
of the free electron and of the electron bound in the H-atom. Half a century 
after its discovery, investigations of Lamb shift still provide a precision test 
of bound-state QED [[| . Now, an experimental uncertainty of 3 ppm in laser 
experiments is essentially smaller then the 10 ppm theoretical inaccuracy due 
to poor knowledge of the proton charge radius. In this situation the study 
of unstable leptonic atoms becomes competitive to the study of hydrogen. 

In a study of the inter-nucleon potential there have been some early expec- 
tations j|] to reveal the off-shell parts of the nucleon-nucleon bremsstrahlung 
amplitude. Fearing and Scherer |5|,|| excluded this possibility by subsuming 
such off-shell amplitudes into redundant terms |7j that can be rotated by 
field redefinitions into so called contact terms. The parity violating anapole 
terms [§ might be one exception, deserving a separate study. 

However, particle physics provides new microscopic interactions leading 
to potentially interesting new contact terms. The most famous one is the 
anomalous 7r°77 coupling. As explained in some detail in || this coupling 
can be viewed as an off-shell effect. On top of the QCD binding of the quark- 
antiquark atom, the triangle quantum loop dominated by far off-shell quarks 
produces an anomalous coupling responsible for the ir° decay. Such manifes- 
tation of the off-shellness motivates us to study the two-photon annihilation 
of atoms in general. 

The simplest "total disintegration" of an "atom" occurs when it consists 
of a particle-antiparticle pair, like in the case of the true QED-atom, positron- 
ium[]. Actually, a comparative study of QED and QCD atoms has been very 
fruitful in the early days of quarkonia. A total disintegration of an atom 
consisting of different fermions is more subtle. It happens on account of the 
flavour changing (FC) processes familiar from weak interactions. While the 
transitions among charged quarks are well known, the lepton-flavour violating 
(LFV) transitions among charged leptons is an open urgent issue, stimulated 
by accumulated indication of the neutrino oscillations. In order to benefit 
from the crossfertilization of different fields we pursue here the comparative 

-"A revival of positronium appeared after the discovery of the QCD atoms, together 
with the recognition that the first proposal of the positronium (termed "electrum" ) has 
been given as early as in 1934 pl| , immediately after the discovery of the positron by 
Anderson. 
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study of annihilation in the heavy-light QED and QCD systems. 

Such a comparative study throws new light on the off-shell nonperturba- 
tive effects of valence quarks, studied first by two of us in the case of the dou- 
ble decays of the K L |[2],[l3|] and B s meson flU}] . Subsequently, this study has 



been continued within the specific bound state models, both for Ki — > 27 [T5|] 
and for — > 27 fTEf . In these papers it was explicitly demonstrated that 



operators that vanish by using the perturbative equations of motion gave 
nonzero contributions for processes involving bound quarks. The purpose of 
the present paper is to elaborate to more detail our more recent study [IT 
which accounts for similar effects for the bound leptons. 



1.1 The relativistic QED atom 

Our starting point is the Lagrangian density for the two fermions in the 
absence of the FC transitions. Thus the light particle (electron e~ of mass 
m) and a heavy positively charged particle (say muon p + of mass M) interact 
only through electromagnetic field, as given by the last term in 

C = C e + - ^F°0F aP - J a A a . (I) 
The Dirac Lagrangian 

A ; d=d -d, (2) 

for a given particle (i = e, fi) leads to the Dirac equations for ipi and ^ 
treated as independent fields 

V><(i + rm) = , (i# - mi)ipi = . (3) 

Imposing the Coulomb (radiation) gauge, V • A = 0, one is able to solve for 
A (eliminate it from the Lagrangian), leading to 

£ = £ e + £, + V ± -B 2 ) + J- A- 1 - f ^r Mr' t) 

2 ; 2 J 4tt \r-r'\ K J 

Here the last two terms can be expressed in terms of p and J components of 
the fermion current 

J a = e(</v7 a Vv - &YVe) • (5) 
The corresponding Hamiltonian, after neglecting the self-energy terms in the 
Coulomb interaction, has the form ]Tj| 

H(x) = H(x) Atom + H (x) Rad + F(x)couiomb-in S t + H(x) int . (6) 



1 a 8 

-7 — m,; 

2 ' dx a 
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where 

#Atom = + H e , (7) 

contains the relativistic fermion contributions 

H eM = J d 3 r H° eM = J d 3 r ^ (x) [ - ia • V + m^]^ (x) . (8) 

The electromagnetic piece splits into the radiation part 

H K&d = l -J &h[El{x)+B\x)} (9) 

containing the relevant electric and magnetic fields 

8 A 

E ± = -—, B = VxA, (10) 
and the instantaneous Coulomb term 

#Coulo mb -inst = ^ / y^J>, t)Je(r'> *) ■ (") 

The relativistic QED atoms can be treated to all orders by solving exactly 
the Dirac equation with a Coulomb interaction. This means solving the Dirac 
equation with V(x) = ^qV^x) (V c denotes the Coulomb potential) 

[i$ + V(x) -rmtyi = . (12) 

Correspondingly, the fermion propagator in external field reads 

[10 + V (x) - m e ]S e F (x, y) = S (4 \x - y) . (13) 

Thus, in contrast to the free-particle propagator, the propagator for the 
bound fermion, 

n,tj n,cr 

should require the sum over all possible excited states, which appear when 
decomposing the fermion field in terms of a complete set of positive and 
negative energy eigenfunctions: 

1> e (x) = £ {b n , a ^(x) + dl a ^_ a (x)} . (15) 
The solutions of the free Hamiltonian 

Hq — Hxtom + -f^Rad + -^Coulomb-inst (16) 
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form a complete set of stationary states \a, N), expressed as a direct product 
of atomic wavefunctions ip a and the photon Fock states 

\a,N) = il> a (r)\N). (17) 

When the interaction is turned on, one should make a replacement 

H ^H = H + H I (t), (18) 

and the pertinent states cease to remain stationary. Their evolution in time 
in practice means that the excited states decay under the influence of the 
QED interaction 

Hj(t) = - J d 3 r [j p (r, t) + J e (r, t)] ■ A(r, t) (19) 

into other states \b,N), where N photons are emitted. 

In addition to the ordinary interaction fll9|), in the next subsection we 
shall consider possible additional interactions Q25D, involving the flavour 
changing fi <-> e transition. This will enable the atom to disintegrate com- 
pletely. The lowest order disintegration requires \i and e overlap, happening 
when / = 0, i.e. a decay from S states. The decay from I > corresponds 
to a cascade down to S state, followed by the decay from there — a higher 
order process which we do not need to consider in what follows. 

1.2 Some motivation for scrutinizing muonium 

There has been a considerable revival of the interest in muonium (Mu = /i + e~ 
system) in view of the very precise measurements in this system. At the same 
time, the theoretical predictions are plagued by the nonperturbative bound- 
state effects. The only known way to achieve the required precision for the 
bound states is by expanding around a nonrelativistic limit. Such methods, 
like Non-Relativistic Quantum Electrodynamics |nj| , Z0 | start from the bound 



state described by a Schrodinger wave function, and build up corrections in 
terms of the relative velocity of the components. 

For the muonium at hand our analysis and results bear a close analogy to 
the correction to the muon lifetime due to muonium formation, reported in 
| |21| |. In this system electron and muon have r.m.s. velocities f3 e = a ~ 1/137, 
and Pn ~ am e /m^ ~ 3.5 • 10 -5 . In terms of these parameters the bound- 
state corrections acquire a form a n (m e /m A1 ) m , where the corrections up to 
n + m = 4 matter in practice. The current world average for the muon 
lifetime measurements Wl 



r M = 2.19703(4) x l(T b s , (20) 



5 



has an uncertainty of only 18 ppm. In order to benefit from an improvement 
the measurement of r M (and thereby of G^) by a factor of 20 (i.e. reducing 
its uncertainty to only ±1 ppm) a knowledge of modification of due to the 
formation of muonium is required. 

The reexamination of the muonium bound state effect [21] showed only 
a tiny effect, the ~ 6 • 10~ 10 correction to the lifetime 

r Mu = T,(l + f?§). 

This negligible overall shift is in contrast to relatively large 0(a—) velocity 



effects on the spectrum |[21|| . In the present work we are pointing out the 
off-shell effects fl53|) in the radiative annihilation of muonium which are in 
between of these two. 

For completeness, let us mention that besides the mentioned radiative an- 
nihilation, there is also the ^-exchange annihilation Mu^ v e v^ (the analog 
of fi~p capture) with the rate 



am, 




r(Mu -> i/ e i/ M ) = 48vr e -\ T(/i+ -> e + v e u„) . (21) 



Still, it leads to a miniscule branching ratio ~ 7 ■ 10 -12 , and moreover is 
restricted to the orthomuonium decay, which is out of our scope here. Let us 
note that some interest in radiative orthomuonium decay might come from 
the three-photon analog decays: the puzzling discrepancy in orthopositron- 
ium (a brief sketch of the recent status can be found in p3j) and the surprising 



suppression in theoretical estimate for Kl — > 37 |[24|| . These three-photon 
decays provide three-party entanglement similar to the one in quantum op- 
tics [[§]. 

Of course, the muonium annihilation involves the LFV transition which 
is a matter of the beyond the standard model (BSM) physics. Such lepton- 
number violating interaction induces simultaneously a fi —>■ e'j transition |26[ . 



so that the unknown details will cancel in the ratio of these two processes. 

For the heavy-light muonium system [i + e~ (where m M = M ^> m e = m), 
the bound-state calculation corresponds to that of the relativistic hydro- 
gen. Thereby we distinguish between the Coulomb field responsible for the 
binding, and the radiation field |27] participating in the flavour- changing 



transition at the pertinent high-energy scale. In this way, the radiative disin- 
tegration of an atom becomes tractable by implementing the two-step treat- 
ment [2j| : "neglecting at first annihilation to compute the binding and then 
neglecting binding to compute annihilation". This factorization of scales 



was introduced for the first time by Wheeler ||29|| . For the muonium atom at 
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hand, the binding problem is analogous to a solved problem of the H-atom. 
In this way we avoid the relativistic bound state problem, which is a difficult 
subject, and we have no intention to contribute to it here. 

The mentioned two-step method is known to work well for disintegration 
(annihilation) of the simplest QED atom, positronium. Generalization of this 
procedure to muonium means that the two-photon decay width of muonium 
is obtained by using 

= |^(0)| 2 |.M( / , + e--, 77 )| 2 

UirMm ' 1 ' 

where |?/>(0)| 2 is the square of the bound-state wave function at the origin. 
After this factorization has been performed, the rest of the problem reduces 
to the evaluation of the scattering-annihilation invariant amplitude M.. In 
the case of positronium this expression will involve equal masses (M=m), 
and the invariant amplitude which for a positronium annihilation at rest has 
the textbook form |30| 



m = ^v s {v2){mh +mh}M P i) . (23) 

Only the antisymmetric piece in the decomposition of the product of three 
gamma matrices above 



{ } 



y^meih - kMO^D* . (24) 



contributes to the spin singlet parapositronium two-photon annihilation. 
This selects (e^ x e^), a CP-odd configuration of the final two-photon state. 
We will see that for muonium annihilation also the CP-even e\ ■ e\ configu- 
ration contributes. 

The paper is organized as follows: In section 2 we consider the quantum 
field treatment of the annihilation process fi + e~ — > 77 in arbitrary external 
field(s). In section 3 we relate the binding forces to the external fields of 
section 2. In section 4 we perform the calculation of B s — > 77 in several 
different QCD models. In addition we consider the related off-shell bound- 
state effects in —>■ K*j decay. In section 5 we present our conclusions. 



2 Flavour-changing operators for fi + e —> 

Augmenting the electroweak theory by LFV enables the one- and two-photon 
radiative decays /i — >• and \i e^yy. Accordingly, the double-radiative 
transition is triggered by two classes of one-particle irreducible diagrams 
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Figure 1: The examples of the one-particle-irreducible diagrams leading to 
the double-radiative flavour-changing transitions. Only the second-row dia- 
grams exist for the leptonic case. 



(Figs. |T]a and b), related by the Ward identities. After integrating out the 
heavy particles in the loops, these one-loop electroweak transitions can be 
combined into an effective Lagrangian fll3| , 



C(e - /i) 7 = B e» vXp F^ (f i D x lp L if,) + h.c. , (25) 

where muon and electron are described by quantum fields \I/ = if)^ and ip = 
ip e . Correspondingly, for B° s — > 27, the involved fields are ip s = s and ipb = b. 

In our case, we do not need to specify the physics behind the lepton- 
flavour- violating transition in (p5]). For instance, the strength B might 
contain Maki-Nakagawa-Sakata |HJ parameters, analogous to the Cabibbo- 
Kobayashi-Maskawa parameters Ackm in the quark sector. 

Keeping in mind that the fermions in the bound states are not on-shell, 
we are not simplifying the result of the electroweak loop calculation by using 
the perturbative equation of motion. Thus the effective Lagrangian fl2"5|) ob- 
tained within perturbation theory splits into the on-shell magnetic transition 
operator £ CT 

£ ff (l7) = B a ^> (Ma ■ FL + ma ■ FR) i/j + h.c. , (26) 
and an off- shell piece Cp |13| 

C F = B F f [(i Jfi — M) a ■ FL + a ■ FR(ip- m)]ip + h.c. , (27) 
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Figure 2: The two-photon contact (seagull) diagram that can be rotated 
away by a field redefinition. 



where a ■ F denotes a^F^, and L — (1 — 7s)/2 and i? = (1 + 7s)/2 denote 
left-hand and right-hand projectors. To lowest order in QED (or QCD) 
Bp = B a = B, but in general they are different due to different anomalous 
dimensions of the operators in (|26|) and (^3). Let us note that the off-shell 
part Cf has zero anomalous dimension [ p^j . 

By decomposing the covariant derivative, \Jj) = i$ — e4, in the off-shell 
operator (^7|), we separate the one-photon piece 

£ F (1 7 ) = B F $[(i 1$ -M) a-FL + a- FR(i$ - m)]i/> + h.c. , (28) 
from the two-photon piece 

£ F (2 7 ) = B F V[-e4tT ■ FL + a ■ FR(-e4)}ip + h.c. . (29) 
The amplitude for the two-photon diagram (Fig. |^) is given by 

A a =iJd 4 xC F (2 1 ) = A L a +A^ } (30) 

in an obvious notation. The single-photon off-shell Lagrangian Cp^j) leads 
to the amplitude with the heavy particle in the propagator 



A h 



\Bp I / d 4 xd 4 yf (y) -ieUv) iS$\y,x) 



[i % —M) a ■ F 1 {x)L + a ■ F 1 {x)R{i^ x - m) 



ip(x) 



(31) 



and a similar amplitude with the light particle in the propagator 



A r 



i Bp //d 4 id 4 t/^(i) 



[i$ x -M)a-F 1 (x)L + 



a-F 1 (x)R(i$ x 



m 



x iS { p\x,y) -ie4 2 (y) ip(y) 



(32) 



The subscripts 1 and 2 distinguish between the two photons. It is understood 
that a term with the 1 <-> 2 subscript interchange should be added in order 
to make our result symmetric in the two photons. 



9 




(a) 



Figure 3: The shaded boxes indicate the combination of the unrotated off- 
shell transition (proportional to Bp) and the on-shell magnetic moment tran- 
sition (proportional to B a ), giving the effective vertex in Eq. ([47]). 



Within the quantum field formalism, the sum of the equations (|30|) , fl3"i"D 
and fl32|) describes the process fi + e~ — > 77, or \i — > 677. 

Let us now be very general, and assume that both particles (e and //) feel 
some kind of external field(s) represented by V( e ) and V( M ), and obey one-body 
Dirac equations 

[i# - V (i )(x) - m(i)]^(i) = , (33) 
for i = e or fi (in general V(j) = 7 a V^\) , and accordingly the particle propa- 
gators Sp satisfy: 

[$ - V {t) (x) - mi }Sf(x, y) = S^(x - y) . (34) 
Our photon fields enter via perturbative QED, switched on by the replace- 



ment — > Dp = <9 M + ieA^ in flTjp. It should be emphasized that A^(x) 
represents the radiation field and does not include binding forces, which will 
in the next section be related to the external fields V(i). 
Now, using relations flT2p and (|IB|) we obtain 

A b = -A L a + AA b , A c = -A% + AA C , (35) 

resulting in a partial cancellation when the amplitudes are summed 

A a + A b + A c = AA b + AA C . (36) 

This shows that the local off-diagonal fermion seagull transition of Fig. § 
cancels, even if the external fermions are off-shell. The left-over quantities 
AA b and AA C involve the integrals over the Coulomb potential and represent 
the net off-shell effect. 

There are also amplitudes Ad and A e which are counterparts of A b and 
A c when Cf(1j) is replaced by C a . The total contribution from our flavour- 
changing Lagrangian (Cp and L a parts) is then given by 

A d + AA b = i //d 4 xd 4 ^(2/)[-ie4 2 (?/) 
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xiSP(y,x)Q(x)^(x), (37) 
represented by Fig. ^a, and a similar one 

A e + AA C = iJJd 4 xd%^{x)Q{x)iS ( F \x,y) 

x -ie4 2 (yj\i{>(y) , (38) 
corresponding to Fig. |b. The operator Q(x) in these expressions reads 
Q(x) = [B a M + B F V M (x)} a ■ F x {x)L 



+ a ■ F 1 (x)R [B a m + B F V (e) {x)} 



(39) 



The result given by Eqs. (|37Q-(|39D can also be understood in terms of the 
following field redefinition. Eq. fll2"D can be obtained from the Lagrangian 



C D (V,i>) =V[ip- V M — M]ty + 4>[ip— V (e} - m]tf> 
Now, by defining new fields 

= $ + B F a ■ F Lip , ijj' = ip + B* F a ■ Fm , 

we obtain 

C D ($, *) + C F = C D (iJ/, + AC B , 



(40) 



(41) 



(42) 



which shows that C F can be transformed away from the perturbative terms, 
but a relic of it, 



AC B = B F #[V M <7 ■ FL + a ■ Fi2V( e )]V> + h.c. 



(43) 



remains in the bound-state dynamics as a contact term. Thus, the off-shell 
effects are non-zero for bound external fermions. Combining ACb and C a , 
we obtain 

ACb + £* = ^Q*P +h.c. , (44) 

where Q is given by (|39|) . This shows how the upper field redefinition rotates 
away the contact term shown in Fig. ^|, leaving us with the result given by 
Eqs. (P)-(H). 



3 Off-shellness in the muonium annihilation 
amplitude 

The preceding section shows how far we can push the problem within quan- 
tum field theory. Up till now we have made no approximations except for 
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standard perturbation theory. Now we apply the obtained results to the 
double radiative annihilation of muonium. Naively, the product $f ip corre- 
sponds to the bound state of /i + and e~, which might be true only for the 
asymptotic free fields. However, relativistic bound state physics is a difficult 
subject, which we circumvent by sticking to the two-step procedure [||[] as 
explained in Section 1. We perform the calculations in the muonium rest 
frame (CM frame of /i + and e~) where we put the external field(s) equal to 
a mutual Coulomb field, — > 70 (where Vc = — e 2 /4iir). In calculating 
the p + e~ — > 77 amplitude in momentum space, we take for V c the average 
over solutions in the Coulomb potential, which is (Vc) = —(ma 2 /2). In this 
way the muonium-decay invariant amplitude acquires the form which is a 
straightforward generalization of the positronium-decay invariant amplitude 



(|23|) in momentum space. 

The amplitudes + AA^ from Eq. (|3"7D, together with A e + Av4 c from 
(|38|), transformed to the momentum space take the form 

M = H^^( Pa ){H^p_p^ 2 + (i < _ > 2)}ti B (p 1 ) J (45) 

where and u e are muon and electron spinors, and e\ 2 are photon po- 
larization vectors. The factor, incorporating the binding in the form of a 
four- vector U a = (p, 0), 

P = {l-x1p)1l&L + x# J £R(l-1p) , (46) 

accounts for the aforementioned factorization of a binding and a decay, and 
is represented by the shaded box of Fig. H: 



M(l - xpj°)a -FxL + ma- F X R{1 - p 7 c 
Here we introduced abbreviations for two small constant parameters, 



(47) 



M ' P ~ mB a ' 1 ' 

in terms of which the sought off-shell effect will be expressed. Note that 



in the effective interaction (|47|), the left-handed part corresponding to V( M ) 
has gotten an extra suppression factor x = m/M in front of the binding 
factor p, in agreement with the expectation that the heavy particle (p + ) is 
approximately free, and the light particle (e~) is approximately the reduced 
particle, in analogy with the H-atom. 

The annihilation amplitude ([45]) can now be evaluated explicitly. The 
usual procedure of squaring the amplitude and using the Casimir trick for 
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converting spinors into Dirac matrices would give us expressions with traces 
of up to twelve Dirac matrices, making the calculation unnecessary extensive. 
It is much easier to proceed by going into the frame in which the muonium 
is at rest and photons are emitted along the z-axis, i. e. 



fci 



/ 




\ 




( 




\ 
































V 


Id 


J 




K - 


-uj 


) 



1 

72 



( °\ 
1 

±i 

V oy 



(49) 



where u = (m + M)/2 w M/2 is the photon energy In this frame fa and ft 
= 1,2) formally anticommute 



^ = c( 7 °± 7 3 )- 7 =( 7 1 ±i 7 2 ) = -^ l 
so we can group them together and calculate 



(50) 



-2u' 



(et, x e*) • fci cr 



* ^-3 



ei a 



V 



€2 ■ e^J; cr 3 — (e£ x e*) ■ k\ e\ ■ e\ — (e* 2 x el) ■ k\ 



a" 



J 

(51) 

It is now easy to multiply this by the appropriate chiral projectors L and 
R, p7° terms, and v^ip^) and u e (pi) spinors. Now, taking into account that 
muonium leading to the two-photon final state is in the spin singlet state, we 
get the result 



M 



-2eB a M 2 



I2M 



m 



[1 — x 2 + xp + x 2 p)e* 2 ■ e\ 



+ i(l + 2x + x 2 + xp - x 2 p)(e* 2 x e*) • fci 



(52) 



In comparison to the expressions and (^) for parapositronium, we notice 
that in addition to e\ x there appears also e\ ■ e*, a CP-even two-photon 
configuration. 

The explicit expression for p depends on some assumptions. As explained 
previously, we use (Vc) = —ma 2 /2 which gives p = a 2 /2 for B a = B F = B, 
which is a good approximation in the leptonic case. 
Eq.(^) finally gives 



2aM 4 



\^)\ 2 \B (T \ 2 {l + 2xp) 



(53) 
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where we have kept only the leading term in p and x. Since the wave function 
at the origin appears as a prefactor, it is not necessary to know the precise 
value of |?/>(0)| 2 ~ (ma) 3 /7T, in order to know the relative off-shell contri- 
bution. Thus, for muonium, the sought off-shell contribution is only a tiny 
correction, 2xp = a 2 m/M ~ 2.6 ■ I0 -7 , to the magnetic moment dominated 
rate. 

We may note in passing that we have checked our results also by the 
direct calculation of the squared Feynman amplitude ( |45| ) on the computer 
using the FeynCalc Mathematica package for algebraic manipulation of ex- 
pressions involving Dirac matrices and spinors Q . Here the explicit 



Lorentz covariance was preserved at all steps of the calculation and the final 
result was in agreement with the one obtained by hand calculation. 



4 Off-shellness in —> 77 

In comparison to a tiny effect in the preceding section, we expect the cor- 
responding off-shellness in a strongly bound QCD system to be significantly 
larger. We also take into account the B F /B a correction in fl53|), when con- 
sidering the B° s — > 77 decay. 

The expressions (^) to fl2"U| ) apply to the b — > 577 induced B° s — > 27 
decay amplitude by simple replacements // — > s and e — > b. Then one has to 
scale the operators Cf,o defined at the Mw scale, down to the S-meson scale. 
The coefficients B F of Cp, and B a of C a , in Eqs. (^) and (|26|), both being 
equal to B at the W scale, may evolve differently down to the /1 = m b scale. 
This difference between B F and B a is due to different anomalous dimensions 
of the respective operators. Within the SM one can write 

B„f = i| A C km Of . (54) 



The coefficient Cj has been studied by various authors p4"|-|58|]. The coeffi- 
cient Cj was considered in fI3| , where at the 6-quark scale we obtained 



|L ~ 4/3 (jx = m b ) . (55) 

Although the off-shell effect for B — > 27 is expected to be suppressed by the 
ratio (binding energy)/mb, it could still be numerically interesting. 

The conventional procedure when evaluating the pseudoscalar meson de- 
cay amplitudes is to express them in terms of the meson decay constants, by 
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using the PCAC relations 



^ l5 b\B° s (P)) 
(0\sy 5 b\B°(P)) 



if B M B . 



(56) 
(57) 



These relations will be useful after reducing our general expression (f45|) con- 
taining the terms with products of up to five Dirac matrices. After some 
calculation we arrive at the expression for the B s meson decay at rest, which 
is analogous to, and in fact confirms our previous relation (^) obtained in a 
different way, 



M 



B 



6 X 



[1 - x 2 + xt + x 2 t) e* ■ e* + 



+ i(l + 2x + x 2 + xt-x 2 t)(e* 2 x e*) ■ k] 



(5f 



Here, the parameter r represents the off-shell effect in the QCD problem at 
hand, and will be more model dependent than its QED counterpart p. With 
the amplitude fl58|) , keeping only the leading terms in r and x, we arrive at 
the total decay width 



a 



M 5 



18m 2 



f 2 \B a \ 2 (l + 2xr) 



(59) 



where by switching off r we reproduce the result of Ref. [39 



4.1 Coulomb-type QCD model 

In order to estimate the value of the off-shell contribution r, in this subsection 
we assume a QED-like QCD model with the Coulombic wave function 
i/j(r) oc exp(— mra c s). Thus we rely again on an exact solution corresponding 
to effective potential V(r) = — 4a e fr/(3r), with effective coupling a e &(r) = 
-(4vr6 ln(rA pot ))- 1 . Here b = (1/8tt 2 )(11 - (2/3) JV/). The mass scale A pot 
appropriate to the heavy-light quark Qq potential is related to the more 
familiar QCD scale parameter, e.g. A pot = 2.23 A^jg (for Nf=3). Within this 
model, we obtain 

r = ■ ( 60 ) 

By matching the meson decay constant fs and the wave function at the 
origin 

N J^El = ffA 2 . |fe(0)|2 = (!^_ (61) 
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we obtain the value for the strong interaction fine structure strength a e g ~1. 
Then, including fl55|) for the QCD case, the correction factor 



xt ss 0.1 , (62) 

is much larger than xp in the corresponding QED case. Correspondingly, 
one expects even more significant off-shell effects in light quark systems, in 



compliance with our previous results [12, 13, 15 



4.2 A constituent quark calculation 

As an alternative to the Coulomb-type QCD model described above, now we 
adopt a variant of the approach in Refs. [[14],|16|]. One might use the PCAC 
relations (|56|)-(|5"T|) together with a kinematical assumption for the s-quark 
momentum, similar to those in Refs. |39|,[|2]]. We assume the bound s and b 
quarks in B° s to be on their respective effective mass-shells. Note that even 



if one is using (|56|) and fl5T|), the amplitude will still explicitly depend on the 
s-quark momentum p s . This is put on the effective mass-shell by using the 
relation = — M s {ki + k 2 )^/M b , where M q = m q + m (for q = b,s) are the 
effective (total) masses, m q are the current masses, and mo the constituent 
mass of order a few hundred MeV. The structure of the amplitude now comes 
out essentially as in (|58| ) with a relative off-shell contribution 

2m n 

xf = « 0.1 , (63) 



of the same order as in fl62|). However, unlike flog), the off-shell effect is now 
only in the CP-odd term (e\ x e^), the square bracket in (|58D being replaced 

by 

e* 2 ■ e* + i(l + 2x + xf)(e* 2 x e*) . (64) 



This may be different in other approaches [|T^] , showing the model dependence 
of the off-shell effect. For instance, potential-QCD models in general, besides 
a vector Coulomb potential, also contain a scalar potential. 



4.3 A bound state quark model 

In our previous accounts [0,|16| we applied a bound state model for B® — > 27. 
Then the potentials Vi in ( |T2"D are replaced by a quark-meson interaction 
Lagrangian 

£*(s, b) = G B bir, s $ + h.c. , (65) 
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B s - 





Figure 4: The two-photon transition amplitude from a contact term (pq) left 
over after the field redefinitions. 



where $ is the B-meson field. In this case, the term Cp can be transformed 
away by means of the field redefinitions: 

s' = s + B F a ■ F Lb , b' = b + B* F a ■ FL s . (66) 

However, its effect reappears in a new bound-state interaction A£$, 

£*(s, b) + C F = C<s,(s', b') + , (67) 

where, after using R-f 5 = R and L75 = — L, 

AU = Bp G B [b'a ■ FL b' — s 1 a ■ FR s']$ + h.c. . (68) 

The two terms in this equation correspond to two contact amplitudes dis- 
played in Fig. |. Also in this case, net off-shell effects are found |T4l,|T6|. 
Further calculations of B — > 27 within bound state models of the type in 
(|5l)D will be presented elsewhere. 

Note that in bound-state models based on heavy-quark effective theory 
the expression ( |5"oD is slightly modified such that the b quark field will be 
replaced by the product of the reduced heavy-quark field and its projector 
P+(v) = (1 +7 • f)/2, where v is the velocity of the heavy quark 



4.4 Link to B d K*-f 

Although we focused till now only to the two-photon processes, the interac- 
tion in (|2"TD contributes to the one-photon couplings as well. 

Actually, we observe that off-diagonal one-photon couplings contained in 
the Lagrangian given by (|39| ) and (|44| ) can be used to calculate the amplitude 
for muonic hydrogen decaying to a photon and ordinary hydrogen, that is, 
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B d \s B * 

d^ivwwww 




K* 

(a) 

Figure 5: Diagrams for Bd — > jK*: (a) The magnetic moment transition 
amplitude, (b) The contact term ( |B"9"D left over after field redefinitions. 

the process pT — > e~ + 7 for both leptons bound to a proton. This is a 
leptonic version of the celebrated 5-meson decay Bd — > K*^. 

As a toy model, one might consider a process "//" — > "e" 7 in an external 
Coulomb field, with and "e" rather close in mass such that the non- 
relativistic descriptions of the "leptons" might be used. The effective — > 
"e" 7 interaction is given in (44). If we assume that (M — m) is of order 



a m, we obtain off-shell effects of order a 2 due to £p, relative to the standard 
magnetic moment term C a . Bigger mass differences give bigger effects, until 
the non-relativistic approximation breaks down. 

Returning to the off-shell bound-state effects in the important Bd — > K*^ 
decay, they can be addressed in the framework of models P3|-fi8| combining 



heavy quark effective theories with the ideas of Nambu-Jona-Lasinio models 
and chiral quark models. 

The ordinary, on-shell transition magnetic moment /^-induced amplitude 
for the Bd — > K*^ is shown in Fig. [|a. 

Now, transforming away the term Hp by the field redefinitions produces 
the new contact term 

= -B F G B s'a ■ F RP + {v) l5 d' <$> (69) 

giving the amplitude displayed in Fig. |S|b. 

The ratio of the off-shell and the on-shell amplitudes in the soft K* limit 
can now be calculated to be 

A(B d -> -g*7)off-shcii _ B F m ^ 

A(B d -> K*7)on-shell B a 2m Q ' 

where m is the constituent mass of order of 200-300 MeV, and rriQ is the 
heavy (b quark) mass. Going away from the soft K* limit, the amplitudes 
will change, but the result (1701) will persist in the leading order. 
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A refined calculation would be desirable in view of the importance of this 
result. An earlier attempt (the preprint version of Ref. reported on 

large off-shell effect in the amplitude which happened to be reduced below 
10 % effect in the rate. 



5 Discussion and conclusion 

The present "atomic" approach enables us to see in a new light the off-shell 
effects studied first for the Kl — > 77 amplitude in the chiral quark model 
12|,|13ll) and subsequently in the bound-state model [[0|. The observation 



that off-shell effects can be clearly isolated from the rest in the heavy-light 



quark atoms |TJJ was still plagued by the uncertainty in the QCD binding 
calculation fllEf . Here, in the Coulomb-type QCD model we are able to 
subsume the effect into an universal binding factor, in the same way as for 
the two-photon decay of muonium in the exactly solvable QED framework. 
It is a quite significant 10 percent effect in the B® —>■ 77 case, whereas in the 
two-photon decay of muonium it is very small (of order 10~ 7 ), but clearly 
identifiable. 

As a byproduct we obtain here also the on-shell amplitude already con- 
sidered in the literature. There is an extensive list of calculations [P^,^2|,pm 
pertinent to the short distance electroweak loop contributions to b — > 577 
which trigger B s — > 77. Comparing our results to the expression (22) of |SS 



we can express Cf in our (p4|) in terms of their coefficient C 

Cr = ^(c+ 2 4), (71) 



aVq V 3 J 

or, numerically, Cf = 0.4 at the B meson scale. Still, there is another 
class of contributions, belonging to the LD regime. For example ]5l]. [j2 



present magnetic moment, 07-type LD effects in B s — > 77 decays in the 
vector meson dominance approach, whereas the other authors |53|-[)I| , though 
with controversial results, estimate the contribution of the charmed-meson 
intermediate states. These seem to be a natural representation for our short- 
distance loops when the loop momenta are below the b quark mass scale. 

Our message is that such small SM effects might obscure possible new 
physics (BSM) signals that are of a comparable size. Without pretending on 
completeness we give some examples that the off-shell effects considered in 
the present paper might hinder possible BSM discovery. 

Let us start with the famous magnetism of the muon, an ideal of the 
precision measurement. At some level the binding effects might become 
relevant. Such (^Bound — 2) effect due to the diagonal one-photon coupling 
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would correspond to the (^Bound — 2) calculated already for a bound electron 



|p6 |. This effect might be interesting in light of a deviation from the Standard 
Model expectation of the order of 10~ 9 recently measured for (g — 2) of the 
positively charged muon [fT^j . Actually, this measurement triggered various 
speculations ascribing this discrepancy to the various BSM effects, the lepton 



compositeness [58 being one possibility. However, there are more direct 



ways to set a bound on the compositeness scale from the flavour-conserving 



processes. For example, there are flavour-diagonal e + e 77 contact terms [59 



^contact = ilPel^DM (j^F^ + ^F^j (72) 

which would lead to a (1 + <5dev) correction factor to the photon angular 
distribution da/dQ in e + e~ collisions. From 

W = s 2 /(2a)(l/At + 1/A*)(1 - cos 2 9) (73) 

LEP200 sets a bound A > 1687 GeV (for A 6 = A 6 = A) at 95 % CL. Thus, 
eventual non-standard BSM physics contribution at LEP energies are highly 
suppressed. 

More promising route to reveal BSM contributions could be provided by 
flavour non-diagonal transitions. Recent evidence for neutrino oscillations 
has renewed interest in charged LFV searches. Among variety of probes 



reviewed in |6(| jj, — > and fi — e conversion (invoking new high energy scale 
M12) seem to be the most promising. Since the effects of new physics are 
expected to enter at one-loop level, these transitions may be parameterized 
by 

c » = 4?ik'»** F ~< (74) 



in order to estimate the sensitivity of the current experimental facilities [pi 

Although the BSM effects might be more pronounced for the flavour- 
changing quark transitions, their discovery might be hindered by the rel- 
atively more pronounced bound-state effects treated in the present paper. 
The off-shell contribution may affect the discovery potential in the radiative 
B meson decays (for example, p2]-|Bl|). In particular, our result (|T0D indi- 
cates hindrance of the BSM discovery potential in the otherwise promising 
B d — ► if *7 decay. 
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